Gauss map and Lyapunov exponents of 
interacting particles in a billiard 



Cesar Manchein and Marcus W. Beims * 

Departamento de Ftsica, Universidade Federal do Parana, 
81531-990 Curitiba, PR, Brazil 



Abstract 

We show that the Lyapunov exponent (LE) of periodic orbits with Lebesgue measure 
zero from the Gauss map can be used to determine the main qualitative behavior 
of the LE of a Hamiltonian system. The Hamiltonian system is a one-dimensional 
box with two particles interacting via a Yukawa potential and does not possess 
Kolmogorov-Arnold-Moser (KAM) curves. In our case the Gauss map is applied 
to the mass ratio (7 = mil 'mi) between particles. Besides the main qualitative 
behavior, some unexpected peaks in the 7 dependence of the mean LE and the 
appearance of 'stickness' in phase space can also be understand via LE from the 
Gauss map. This shows a nice example of the relation between the "instability" of 
the continued fraction representation of a number with the stability of non-periodic 
curves (no KAM curves) from the physical model. Our results also confirm the 
intuition that pseudo-integrable systems with more complicated invariant surfaces 
of the flow (higher genus) should be more unstable under perturbation. 
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1 Introduction 



The Gauss map [Tf2f3] is a chaotic map which generates the integers from 
a simple Continued Fraction (CF) representation of a real number. CFs have 
been used in several different scientific context like, for example, the renormal- 
ization group theory [lf5|6] . expansion technique applied to a model for Bloch 
electrons in a magnetic field |7J, stability of elementary particles [H] and their 
mass ratio representation 09]. In the context of nonlinear dynamics it has 
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been applied to compute stable and unstable directions of maps [ID] , to ap- 
proximate irrational winding numbers for the KAM tori [11] and to determine 
critical parameter values at which the KAM tori break [T2ll3|ll4] . 



In this paper we discuss the stability of perturbed tori in a different situation. 
The unperturbed system is the problem of two particles in a ID-Box inter- 
acting via Hard Point-like Collisions (HPC). Depending on the mass ratio 
between particles, the system can be integrable (Invariant Surface (IS) of the 
flow is a torus with genus g = 1), pseudo- integrable (IS has a more complicated 
topology with g > 1) and ergodic (see [151161117] ). Since the HPC case is lin- 
early unstable, the LE is zero for each mass ratio and little is known about the 
stability of the IS. Since the Hamiltonian is not differentiable at the boundary 
and at the collisions, no KAM curves exist for any perturbation. Classical [15] 
(quantum [TBlfTB] ) results strongly suggest that the classical dynamics (level 
statistics) is more unstable (Gaussian Orthogonal Ensemble (GOE) distribu- 
tion) for the pseudo-integrable than the integrable cases. Here the perturbed 
system is obtained by changing the HPC to a Yukawa interaction (YI). Recent 
classical results [18] also suggest that pseudo-integrable IS from the HPC are 
more unstable than the integrable ones when the YI is turned on. We show 
here that the chaotic dynamics in the perturbed system, the ID-Box with YI, 
is directly related with the chaotic property of the number related to the mass 
ratio 7 = mj/mi between particles. To do this we compare the LE from the 
infinite continued fraction representation of 7 with the maximal LE from the 
two interacting particles in the ID box. In fact, we show that the qualitative 
behavior of the LE in the YI case can be reproduced qualitatively using the 
LE from the Lebesgue Measure (LM) zero Periodic Orbits from Gauss Map 
(POGM). We also show that the dynamics in the Yukawa case is more unstable 
when POGM are closer to non-POGM. In these cases 'stickiness' (tendency 
of orbits to get trapped) appears more often in phase space. This is shown by 
using the most probable Lyapunov exponent, proposed recently [TS] as a very 
sensitive tool to probe globally details in phase space dynamics. 

It is well known from the KAM theorem that in two-dimensional systems the 
torus surface does not exist if a = uji/uj 2 lies in a region for rational numbers, 
where W\ , w 2 are the frequencies of the unperturbed problem. The region for 
rational numbers increases as the perturbation parameter increases. There- 
fore, only those irrationals that are hardest to approximate by rationals will 
yield the KAM surface. The residue criterion [13] establishes a correspondence 
between the existence of a KAM curve and the stability of the periodic or- 
bits that approximate it. In other words, the stability of periodic orbits very 
close to non-periodic orbits (KAM curves), allows to make some statements 
about the destruction of the KAM curve. In our case the situation is different, 
the LE calculated for the periodic orbits with LM zero from the Gauss map, 
i. .e. the non-periodic curves from the physical model (not KAM curves), al- 
lows us to get some insight about the stability of the IS from a linear unstable 
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system with LE zero. Therefore we have a direct relation between the invariant 
surfaces with zero LEs and their stability under perturbation. 

In section [2] the main results from the ID-Box HPC are summarized. Section [3] 
analyzes those properties of the Gauss map which are relevant for the purpose 
of this paper. We calculate explicitally the LE for the POGM with LM zero 
which differ from the ergodic result. Section H] introduces the smooth Yukawa 
interaction between particles and the LE distribution is calculated. The re- 
lation between the mean LE from the physical model and the LE from the 
POGM is demonstrated. The paper ends with conclusions in section [5j 



2 Two particles in a ID-Box with hard point-like collisions 



Two particles in a ID-box with HPC can be treated as a particular case of the 
motion of three particles on a finite ring [T9|20] . which can be mapped onto 
the motion of a particle in a triangle billiard [2T] . The whole dynamics can be 
monitorated by changing the angles of the triangle billiard. These angles are 
functions of the masses ratio between particles. It is possible to show [22J that 
the dynamics is non-ergodic if 9 is a rational multiple of tt, where 



l-m 2 /m 1 1-7 

cos (0) = — — = — — . (1) 

l + m 2 /m 1 1 + 7 

Writing 9 = —tt, where m and n are integers, at most 4n distinct velocity 
values occur. These are the periodic orbits from the problem. Although there 
are infinite mass ratios which give rational values of 9/tt, some of them are 
special: the integrable cases [23] 7 = 1,3 (or 1/3), which have 9 = ~7r and 
9 = |7r (or 7r/3), respectively. These are the cases for which the genus g = 1 [21] 
(the IS of the billiard flow is a torus). For all other rational 9/tt the dynamics 
is pseudo- integrable [2U], the invariant flow is not a torus and gets more and 
more complicated as g increases (1 < g < 00) [17] . It was shown [TB] that for 
small g the spectral statistics is close to semi-Poisson and it approaches the 
GOE statistics when g is increased. Table [1] shows some rational angles for 
the right triangle billiards with g < 3 and their relation with the mass ratio 
7 and 9/tt obtained from Eq. (TjQ). For genus g = 2, for example, the values of 
9 are: jTT, ^tt and the mass ratio are respectively 7 ~ 0.17 and 7 ~ 0.11. On 
the other hand, when 9 is an irrational multiple of tt, the velocities become 
uniformly dense |25j in velocity space. As a consequence, it is at least possible 
for the two-particle with HPC to be ergodic in velocity space. 

Since the ID-Box HPC is linear unstable [26], all LEs are zero and little is 
known about the stability of the above mentioned ISs. In this work we propose 
to analyse the stability of the ID-Box HPC using results from the Gauss map. 
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genus (g) 


p/q 


6/tt 


7 


l/ 7 




x km (CF) 


1 

Integrable 


1/4 


1/2 


1 


1 


7.37 


[2, 800 


-1/2 


1/6 


1/3 


1/3 


3 


7.78 


[3, 800 


~ 1/3 


2 

Pseudo-Int. 


1/8 


1/4 


0.17 


5.89 


8.07 


[4, 800 


~ 1/4 


1/10 


1/5 


0.11 


9.47 


8.30 


[5, 800 


- 1/5 


3 

Pseudo-Int. 


1/12 


1/6 


0.07 


13.9 


8.47 


[6, 800 


- 1/6 


1/14 


1/7 


0.05 


19.2 


8.64 


[7, 800 


- 1/7 



Table 1 



Some rational angles p/q (see [16]) for genus g < 3 from the triangle billiards and 
their relation with the mass ratio 7 and 6/tt obtained from Eq. ([1]). A (x/^soo) is 
the LE obtained from period-2 orbits of the Gauss map calculated very close (on 
the left) to 6 /n. The last row shows the CF representation of the Xk soo used to 
calculate the LE. 

We will show that, associating the irrational values of 6/tt (and therefore 7) 
with properties of the LE for the LM zero POGM, we are able to make some 
statements about the stability of the ISs with LE zero from the HPC case. 



3 The Gauss map 



In this section the main properties of the Gauss map will be reviewed and 
some numerical calculations for the LE will be performed. For more details 
we refer to the works of Corless et al [2|3] and references therein. The Gauss 
map for x in the interval (0, 1) is given by 



G{x) 



1 

x 



r 

.X _ 



(2) 



The notation [.] means to take the fractional part. The LE exponent can be 
calculated from 



X G (x) = lim -In ( TT \G'(xi)\ \ , (3) 

whenever this limit exists, where G'{xj) = dG(xj) / dx^. The Gauss map gen- 
erates the numbers n\, 712, n$, . . . from the simple CF representation of a real 
number 
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where the are all positive integers, except no which may be zero or neg- 
ative. Here we consider n Q = 0. The CF is represented in the form x = 
no + [ni, n 2 , n,3, . . .]. For rational values of x the sequence of n, is finite and 
the LE from Eq. [3] cannot be calculated. For irrational x the sequence is in- 
finite. For irrational quadratic values of x, the sequence of rii are periodic, as 
exemplified below. The Gauss map is ergodic and for almost all initial condi- 
tions the LE may be calculated explicitly by X G (x) = 7r 2 /61og2 = 2.3731 . . .. 
Periodic orbits in the Gauss map occur when a sequence of integers repeat. A 
fixed point has the property x ni = [n 1; n 1; n 1; . . .], which is written in a more 
simplified manner as [n"i]. For period-1 periodic orbits the LE can also be cal- 
culated explicitly, however it is different from the almost-everywhere values 
(in the Lebesgue sense). It can be calculated from A G (x ni ) = 2ln(l/x ni ). For 
the golden- mean number x = — l)/2, for example, we have n\ = 1 and 
X G (x\) = 2ln(l/x ni ) = 0.96 . . . which is the lowest LE. All initial conditions in 
(0, 1) which have the property x = [n 1; n 2 , n 3 , . . . , n/., 1, 1, 1, , . . .], have the LE 
equal to ~ 0.96. Period-2 orbits have the form x = [n%, n 2 , n 1; n 2 ■ ■ .] = [ni, n 2 ], 
period-3 orbits have x = [ni, n 2 , n^}, and so on. Note that for any rational x 
the limit j — > oo from Eq.Q does not exist and the LE cannot be determined. 

Now we give a numerical summary of the above results for the LE of period- 
1 and period-2 orbits in the Gauss map. If we iterate the Gauss map many 
times (j — > oo) using arbitrary initial conditions between the interval (0,1), 
all orbits have LEs equal to ~ 2.3731, which is calculated using Eq. (J3j). This 
is the ergodic result and it is shown by the solid line in Fig. [TJ Note that in this 
numerical procedure, due to numerical errors it is impossible to get exactly 
the POGM when j — > oo. As a consequence the LE is always ~ 2.3731. 
However, for low period-p orbits we can calculate the LE which differs from 
this ergodic results, since just few iterations are needed and numerical errors 
do not have time to propagate. We use the following procedure, we chose a 
period-p POGM map, determine the corresponding initial value of x and then 
calculate numerically the corresponding LE A G (x) using j = p. For period-1 
orbits we have for n\ = 1, 2, 3, 4: 





= 0.618033989.. 


. = [1,1,1,.. 


•] = m, 


\ G (xi) ~ 0.96, 


x 2 


= 0.414213562. . 


. = [2,2,2,.. 


•] = [2], 


A G (x 2 ) ~ 1.76, 


x 3 


= 0.302775638 . . 


. = [3,3,3,.. 


•] = [3], 


A G (x 3 ) ~ 2.38, 


x^ 


= 0.236067977.. 


. = [4,4,4,.. 


■] = % 


A G (x 4 ) ~ 2.88. 



14.0 




Fig. 1. Lyapunov exponents for period-1 and period-2 orbits from the Gauss map 
in the interval (0, 1). 

We observe that using these value of Xk {k = 1,2,3,4), after one iteration of 
the Gauss map the LEs differ from the ergodic result. The LEs X G (xk) for 
period-1 points are plotted as circles in Fig. [TJ They are exactly on the dashed 
line, which is the curve lnG'(x) for just one iteration. Points along the dashed 
line are the LE A G (a;) from the Gauss map only when x = x±, X2, ■ ■ ■■ Later we 
will explain the reason why we plotted the dashed line. 

We can also calculate the LE for period-2 orbits. Some explicit example are 



%2,1 


= 0.366025 . . 


•~12,1], 


A G 




~ 1.32, 


%3,l 


= 0.263762 . . 


■ -[ML 


A G 




~ 1.56, 


^4,1 


= 0.207106.. 


■~K 1], 


X G 




~ 1.76. 



We generated the values of x = [n 1 ,n2] for Period-2 orbits using all com- 
binations of ni,n 2 . These points are marked as crossed points in Fig. (TJ 
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where some corresponding CF representations are shown. Different sequences 
(branches) of x values are observed. Sequences start (from below) at [T], [2, 1], 
[3, 1J, . . . [k, 1] ... , [800, 1J. We could increase the last value of k but it is not 
relevant for the discussion here. In all simulations the last k used was k = 800. 
In Fig.[T]we showed only the 15 first sequences, they come closer and closer as k 
increases. Inside each sequence the CF representation changes its second num- 
ber. For example, the first sequence, on the right, starts at the Golden mean 
[T] with the lowest LE for this sequence and the LE increases with period-2 



orbits [1, 2], [1, 3J until [1, 800J, which has a LE ~ 6.69. The next sequence 



of period-2 orbits is [271], [272], [273], . . . , [27800]. The last point from this 
sequence has a LE equal ~ 7.37. In fact, for X2,k=oo = [2, ooj, which is exactly 
equal 1/2, the limit of ([3]) does not exist and the LE cannot be calculated. 
Using the periodic orbit £2,800 = [2, 800] ~ 1/2 we are allowed to calculate 
the LE A (x2,8oo) ~ 7.37 very close (on the left ) to the non-periodic orbit 
x = 1/2. 

The next sequence (see Fig. [[]) starts at [3, 1] and ends at x = [3,800] ~ 1/3 
with LE ~ 7.78. The subsequent sequences converge to [4, 800] ~ 1/4 with 
A G (x 4 ,8oo) ~ 8.07, [57800] ~ 1/5 with A G (x 5 ,80o) ~ 8.30 . . ., [800, 800] ~ 0.0012 
with A G (x8oo,80o) ~ 13.4. It is important to observe that the LE increases more 
and more as k increases. For all these points the LE calculated for the POGM 
differ from the ergodic result ~ 2.3731. While such POGM have LM zero and 
may not be relevant for the Gauss map itself, we will show they contribution 
in a physical problem. We just need to relate the Gauss map and the ID-Box 
HPC from section [2] through x = 9 /7r. 



4 Yukawa interaction, results and discussion 



In order to study the stability of the ISs from the HPC case, we need to apply 
a perturbation on the system. Therefore, we assume now that the interaction 
between particles is given by the Yukawa potential 



6 _ 

V{r) = V , (5) 

r 

which has strength Vq and the parameter a > gives the interaction range 
ro = 1/a. The classical dynamics of this problem was already analyse for equal 
masses j27j and for mass ratios 7 in the interval (1,4) [T8] . 

Using the above interaction, we calculated the finite-time maximal LE as a 
function of the mass ratio in the interval (0, 4). Results are shown in Fig. [2]^a) 
for the distribution P(A t , 7) of the finite-time largest Lyapunov exponents [28] 
A t and Fig. E](b) for the mean (At) (see solid line). The mean LE decreases 
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from roughly 1.18 for 7 ~ 0.0 to 0.59 for 7 = 4.0. An unexpected pronounced 
peak is observed at 7 ~ 1.0. When the value of 7 decreases, the mean LE 
exponent presents a minimum at 7 ~ 0.85 and then increases again until 
another unexpected lower peak close to 7 ~ 0.11 [better seen in Fig. [2fe,)]. The 
question now is: what is the origin of such peaks? Why are there special values 
of 7 where the motion is more chaotic? 




0.0 P" ■ 1 ■ 1 1 ■ h 

0.0 1.0 2.0 3.0 4.0 

Y 



Fig. 2. a)Finite-time distribution of the largest Lyapunov exponent P(Aj, 7) calcu- 
lated over 400 trajectories up to time t = 10 4 and for tq — ► 00. With increasing 
P(A(,7) the color changes from light to dark (white over yellow and blue to black) 
and b) mean Lyapunov exponent (solid black line), normalized distribution -Pa (7) 
of the most probable Lyapunov exponent (solid gray line) and results for the LE 
from the Gauss map from Fig. [TJ 
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Figure [2(b) also shows the results for the 7 dependence of the LE from the 
Gauss map showed in Fig. [TJ We used relations 6 = ttx and j(x) = ^ 
from Eq. ([T]) to calculate the correspondence between 7 and x from the Gauss 
map. Note that in this interval of 7 only the Golden mean point [T] remains 
from the first sequence (on the right) from Fig. [TJ All other points from this 
sequence have 7 > 4.0. First observation is that the main qualitative behavior 
of the LE from our physical model follows the curve calculated for the first 
iteration from the Gauss map (see dashed line). In other words, the main 
qualitative behavior of the mean LE from the ID-box with YI follows the 
period-1 POGM. Second observation is that the pronounced peak observed at 
7 ~ 1.0 is very close to the point 7(0; = [2, 800J ~ 1/2) which is a period- 
2 orbit from the Gauss map. We have to remember that this point is very 
close to the non-periodic orbit x — 1/2 from the Gauss map where the LE 
cannot be calculated. The other pronounced (lower) peak at 7 ~ 0.17 is close 
to the period-2 orbit 7(0; = [4, 800] ~ 1/4) from the Gauss map. As the value 
of 7 decreases, the LE from the ID-box with YI follows all peaks obtained 
from the POGM 2(7) = [M00] (k = 2, 3, ... , 800). The important point is, 
the LE from our physical model follows the behavior of the LE from the LM 
zero POGM. Therefore, both pronounced peaks at 7 ~ 1.0, 0.17 are probably 
related to signatures from the periodic orbits of the HPC, for which the LE 
cannot be estimated using the Gauss map. 

Another interesting feature appears if we calculate the change of the width of 
P(A t ,7) around the most probable Af defined through 



<9P(A t , 7 ) 
dA t 



= 0. (6) 

A t =A? 



This quantity, called -Pa (7), has been proposed [18] as a sensitive measure of 
'stickness' in phase-space which are a consequence of the existence of regular 
islands. Each time this quantity has a minimum, the mixed phase-space of 
a system is expected to have more trapped trajectories. Otherwise it has an 
"ergodic-like" motion, then mostly all initial conditions converge to the same 
finite time LE. For more details about this quantity, some examples and moti- 
vations, we refer the readers to [118] • This quantity is plotted in Fig. Mjo) (see 
gray curve with strong variations). Clearly three minima are observed close 
to 7 ~ 0.25,1.0,3.0 where the dynamics in phase-space has more trapped 
('sticky') trajectories. As shown in another work [18], the values 7 ~ 1.0,3.0 
can be related to the integrable cases (genus g — 1, see Table [lj from the HPC 
case. Here we show additional results for mass ratio in the interval (0.0, 1.0). 
Besides for 7 = 3.0, it is very interesting to observe that -Pa (7) has a min- 
imum for all points for which 7(2 = [fc,800]) (k = 2, 3,..., 800) from the 
Gauss map. As a consequence, -Pa (7) has a minimum in the extended interval 
7 ~ (0.0,0.4), where these points come closer and closer. Since the points 
7(x = [k, 800]) (k = 2, 3, . . . , 800) are very close to the non-POGM (which are 
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the PO from the model), it suggest that the minimum of -Pa (7) are due the 
trapped trajectories reminiscent from the PO from the physical model (with 
HPC or YI). 

Above results allow us to make some statements about the stability of the 
linear unstable dynamics from the HPC case, i. e., the stability from the 
different topological surfaces present in the HPC case. Table [1] gives some 
examples of the relation between the values of x = 6/ir with the LE \ G (xk,8oo) 
from the Gauss map calculated very close (on the left) to these points. For 
example, 6/71 = 1/2(7 = 1) is an integrable case from the HPC problem with 
genus g — 1 (torus). The LE for the POGM calculated at 7(2 = [2, 800]) gives 
A G (x2,8oo) ~ 7.37. We are proposing that this value of the LE gives a possible 
"degree of instability" of the torus with genus g = 1. We observe in Table 
[1] that as we increase g, the corresponding values of \ G (xk,soo) also increase 
and 7 decreases. Therefore, we expect that for higher values of g the invariant 
surfaces from the HPC case are more unstable. This is verified for the problem 
considered in this paper (see Fig. [2]), where the mean LE increases as the genus 
g increases following values from Table [H 



5 Conclusions 

While the KAM theorem [TT] makes some statements about the existence 
of non-periodic orbits (KAM curves), the residue criterion [T3] establishes a 
correspondence between the existence of a KAM curve and the stability of the 
periodic orbits that approximate it. Such criterion can be used to determine 
for which parameter of the model the KAM curves may be destroyed. Here 
the stability of periodic orbits with LM zero from the Gauss map (i. e., the 
non-periodic curves from the physical model) allows us to get some insight 
about the stability of ISs (no KAM curves) from a linear unstable system 
with LE zero . The linear unstable system considered here is the ID-box 
with two particles interacting via HPC. The stability of the IS is obtained by 
calculating the LE (via Gauss map) from the CF representation of the masses 
ratio. By perturbing the IS with a Yukawa interaction between particles, we 
observe that the LEs follow qualitatively the LEs from the Gauss map [see 
Fig. [2b)]. Only periodic orbits with LM zero from the Gauss map seems to 
be relevant. Additionally, the two more pronounced peaks in the mean LE 
[see Fig. |2^)-b)] at 7 ~ 1.0 and ~ 0.17 are explained with results from the 
Gauss map. They are probably related to signatures from the periodic orbits 
of the HPC case, for which the LE cannot be estimated via Gauss map. This 
shows a nice example of the relation between the "instability" of a simple CF 
representation of a number to the LEs from a physical model. We also show 
that pseudo-integrable systems, where the IS has a higher genus g, are more 
unstable under perturbations. This is easy to see from Table [T] and Fig. [2b). 
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As the value of g increases, the corresponding value of 7 decreases and the LE 
increases. Moreless, we were able to show that 'stickness' effects are present 
each time the mass ratio 7 is close to POGM and therefore, close to the non- 
periodic orbits from the ID box with HPC This was quantified by calculating 
Pa (7), which is a sensitive measure of trapped trajectories in phase-space. 
Each time -Pa (7) has a minimum [see gray line in Fig. [2b)], trapped trajectories 
in phase-space are expected . 
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